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In the framework of mean-field approximation the dynamics of Bose- 
Einstcin condensates can be described by the hydrodynamic-like equations. 
These equations are analyzed here with account of mutual interaction be- 
tween condensate and non-condensate atoms. The Lagrange invariants and 
freezing-in invariants of such a system have been found. This allows to get 
some necessary conditions for creation of an atom laser with controlled pa- 
rameters of the beam. Particularly, the atom laser beam can carry quite 
well defined angular momentum. This can be practically realized in the 
most simple case, when the vorticity of condensate appears to be a freezing- 
in field. The optimal conditions for a source mode regime are found out in 
the paper. 

PACS numbers: 03.75.Fi 05.30.Jp 47.37.-Fq 

1. Introduction 

One of the youngest, but very rapidly developing field of modern low 
temperature physics is the problem of trapped Bose gases at very low tem- 
peratures, when Bose-Einstein condensation occurs. For recent reviews of 
this topic in general see Refs. [1, 2, 3, 4]. In last few years a consider- 
able progress has been achieved in understanding of the dynamics of Bose- 
Einstein condensates (BEG) both at T = [1] and at finite temperatures 
[2, 5, 10]. However, there many questions remain in the latter case yet 
(see discussion in [11]). To realize the full potential of recent developments 
in BEG physics, and to analyze adequately the data of the experiments in 
a quantitatively meaningful manner, a detailed understanding of the BEG 
dynamics in all its aspects is required. 

A detailed investigation of the dynamics (and, particularly, collective 
processes) of BEG at various temperatures is extremely important in view 
of the possible applications of unique features of such systems. In all like- 
lihood the most intriguing experimental project associated with trapped 
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atomic gases is the so-called 'atom laser' or, in other words, highly coherent 
atomic beam generator [12, 13] (for a review see also [4]). The proposed ex- 
perimental configurations should satisfy a number of basic criteria in order 
to be called an atom laser. Of course, the high phase coherence of atomic 
beam is required first of all. However, the highly dispersive nature of the 
BEC suggests that the spatial focusing and stability of cross section of a 
BEG beam will present possibly more of a problem than encountered in the 
process of focusing laser light [14]. Very recently the quasi-continuous atom 
laser has been constructed [15]. Furthermore the same group demonstrated 
the successful atom optical manipulation, such as reflection, focusing and 
the beam storage in a resonator [16] (however, see the earlier successful 
works [41, 42] based on rather different principle). Thus, another interest- 
ing task appears, namely how to create the atom laser with well controlled 
characteristics. Solution of such a problem can bring atom laser closer to be 
a useful tool in various potential applications. The present paper solves one 
of the problems in this direction. Particularly, we consider the possibility 
for atom laser beam to carry the predesigned angular momentum. 

The paper is organized as follows. In section 2 I consider the quasi- 
hydrodynamic approach to a trapped Bose-gas below the BEC transition 
temperature. Section 3 is devoted to the description of the method for con- 
struction of Lagrange invariants and freezing-in fields by means of the gauge 
transformation of the BEC quasi-hydrodynamic equations. The applica- 
tions of such invariants to the creation of atom laser beam with predesigned 
parameters are considered in Section 4. 



2. Quasi-hydrodynamic equations for the condensate atoms 

The starting point in description of BEC is the usual Heisenberg equa- 
tion of motion for the quantum field operator ip{r,t) 



at 



i>{r,t),H] = ^-^A + U('^-P\r)\^P{r,t)+g\i>{r,t)\^{r,t), 



(1) 

where [/'*''"*'^(r) is the confining potential, the explicit form of which is 
not essential for us here. In equation (1) we assumed s-wave short-range 
interatomic interaction with a strength g = Airah^/m (a is the effective 
scattering length). As usually, set 

i;{r,t)=^{r,t)+iP{r,t), (2) 

where $(r, t) = (■0(r, i)), ip(r,t) is the non-condensate field operator. Tak- 
ing an average of equation (1) with respect to a broken symmetry nonequi- 
librium ensemble, we come to the equation for the condensate wave function 
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[10] 

dt I 2m 



gm{r, t)$*(r, t) + 5(^+(r, t)V5(r, t)^(r, t)), (3) 
where nc{r,t) = |$(r, i)p is the local density of atoms in the condensate, 

n(r,i) = (V'+(r,#(r,t)), (4) 

is the non-equilibrium non-condcnsatc density. Equation (4) involves also 
the anomal non-condensate density rn(r,t) = (^'(r, t)';/;(r. f)) and the three- 
field correlation function {\ip{r,t)\'^ip{r,t)) . The appearance of two last 
terms in (3) is a consequence of Bose broken symmetry in the system. 

The earlier approaches to the equation (3) were based on the assumption 
that all atoms are in the condensate. In this case the so-called Gross- 
Pitaevskii [17] equation appears: 

'^^^dT^ = |-|^A + C/(*™^)(r) +5ne(r,t)| $(r,t), (5) 

which can be conveniently rewritten in terms of the local condensate density 

nc{r,t) = Mr,t)f (6a) 

and local velocity 

Ve(r,t) = -V0(r,t). (66) 

m 

Here ^(r, t) is the phase of the condensate wave function 



$(r,t) = ytM)exp(i^(r,t)). (7) 

It should be noted that the analogous nonlinear Shrodinger equation can 

be obtained quite rigorously in the case of high density, but weak enough 
point interaction [18]. However, in that case the physical sense of is not 
clear enough. So, using (6) and (7), we can present equation (5) as the set 
of two following equations 

^ + V(ncVc) = 0, (8a) 
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where 



2m -s/ric 



It is remarkable that the equations (8) are hydrodynamic looking. This 

fact somewhat confusing on the face of it. Really, equations similar to (8) 
was obtained phenomenologycally to describe a superfluid component of 
liquid helium [19, 20] which is a strongly interacting many-particle system. 
Nevertheless we should keep in mind that in Bose-Einstein condensed state 
we deal with the single condensate wave function <I>(r, f), which allows a 
strong analogy with the order parameter in superfluids.^ 

The next step is to extend the preceding analysis to finite temperatures 
where there is a large fraction of atoms outside of the condensate. In this 
case we need two equations to be used. While the condensate wave function 
can be described as earlier by the Gross-Pitaevskii equation, the distribution 
function of the non-condensate atoms obeys a kinetic equation, which must 
take into account the collisions of both types: non-condensate atoms with 
each other and their interaction with a condensate. The quite rigorous 
derivation of the corresponding collision integrals can be found in Appendix 
A of paper [10]. After some mathematics we come to the corresponding 
hydrodynamic-like equations 

^ + V(ncVc) = r(r,t), (9a) 



where 



r(r,i) = - /(^-^[/(P'^'*)]' (10) 



J [/(p,r, t)] is the collision integral corresponding to the collisions between 
condensate and non-condensate atoms, which functionally depends on the 
distribution function /(p,r, t) of excited atoms. So, function r(r, t) is the 
characteristic rate of the atoms exchange between condensate and non- 
condensate. New chemical potential fi in Eq. (9b) is now defined by the 
relation 

^ = + t/(*™?')(r) -|-5nc(r,i) + 2gnir,t), (11) 



^ Note, there is no complete microscopic theory of supcrfiuids until now. We only know 
that the dynamics of superfluid component in superfluids can be well described by 
the hydrodynamic equations similar to (8) . 
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where 

is the density of non-condensate atoms. Equations (9) describe the dynamics 
of BEC hke an 'effective fluid' with varying density. Term r(r,t) play the 
role of inhomogcneous and nonstationary source. 



3. Hydrodynamic invariants 

As we reahzed in the previous section the evolution of BEC in the frame 
of reasonable approximations can be described by the equations, which are 
similar to the hydrodynamic ones. Classical equations of ideal liquid have 
quite a number of invariants. Except of ordinary integral invariants there are 
local invariants as well. The Lagrange invariants and freezing-in invariants 
are most important ones. Lagrange invariants are conserved along the 'liquid 
particles' trajectories, while the freezing-in invariants are used in reference to 
the fields frozen into a liquid, i.e. the corresponding physical quantity (field) 
vanishes in a frame which moves with the fluid. In the papers [21, 22, 23] 
a wide class of invariants was found. Furthermore, the authors of Refs. 
[21, 22] proposed the method of obtaining new invariants on the basis of 
already known. The recent paper [24] was devoted to construction of the 
invariants of super fluid hydrodynamic equations by means of their gauge 
transformation [25] . This method is very attractive because after the gauge 
transformation the presence of many additional invariants becomes obvious. 

The idea of gauge transformation [24, 25] can be modified to be helpful 
in our case, i.e. BEC dynamics. Really, the local condensate velocity is 
defined by the equation (6). If we wish the condensate wave function to be 
single-valued, then the bypassing along vortex line must lead to the change 
of a phase by the integer of 27r. To be so, we should do a cut. If the leap 
of phase on the bank of cut is proportional to some new function, say Uc, 
then we can do the following gauge transformation 

Vc = -V9 + Va + Uc, (12) 

where a is some gauge function. Gauge of the fields should be done by the 
equation for — Va and by the initial conditions. After the substitution 
of equation (12) into (9b) we get 



{I + (VeV)}.. = + |(a-^)}-(VeV) 



da _ 39 
dxi dxi 
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It can be easily tested that if the gauge function obeys the following equa- 
tion: 

+ ^^^^ 
then equation (13) becomes 

{| + (VeV)}.. = -.e.^. (15) 

The gauge of field Uc is determined by equation (14) and by the initial condi- 
tion for 9 or Uc- The scalar product of the field and the flux line element 
61 behaves like a mass element, which is conserved along any trajectory [24]. 
Direct test show that 

{^ + (VcV)|(Ue<51)=0. (16) 

So, the quantity UcSl is the Lagrange invariant. It should be noted 
that in this case the vorticity w = rotvc becomes freezing-in field [26]. 
Using the analogy with superfluid hydrodynamics wc believe that in BEC 
rotVc = everywhere except the axes of vortices. So, we come to the 
following conclusion: if there were vortices in the BEC initially and the 
invariant Uc6\ is conserved, then there will be the given conserved vorticity 
in any frame moving with the flux lines in future. 

The possible existence of vortices in BEC has been under extensive dis- 
cussion for a rather long time (see e.g. [28, 29, 30, 31, 32, 33, 34]). And 
finally they were recently obtained in the experiments [35, 36, 37, 38, 39]. 

In the next section we consider the consequences of this conclusion for 
possible experimental realization of an atom laser. 

4. The stability of atom laser beam 

A number of atom laser schemes have been proposed during last few 
years. Evident progress is already achieved in the realization of pulsed 
lasers using a matter-wave splitter based on radio frequency (rf) transitions 
[15, 40] and optical Raman transitions [41, 42]. Such schemes, however, 
have several shortcomings, the main of which is the difficulties in achieving a 
continuous refilling. Another schemes, which allow to create the continuous 
wave atom laser can be clearly divided into two distinct classes: optical 
cooling [43, 44, 45] and evaporative cooling [46, 47, 48, 49]. Both models 
are based in principle on the same idea: some source supplies atoms to an 
upper-lying mode of an atom trap. This source mode is coupled to the 
ground state mode (condensate) via a particular cooling mechanism. It is 
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hoped that the macroscopic population in this ground state mode, or laser 
mode, can be built up and coupled to outside world to produce the laser 
output. Independently on particular model, cooling process or, in other 
words, process of increase of the atom population in the ground state mode 
must satisfy the main criterion: the uncontrolled perturbation of ground 
state atoms should be minimal. 

I shall not keep myself in the frame of the particular experiments (even 
successful, such as [41, 42] or [15, 16]). Below both the possible situations 
are considered. 



4.1. Continuous m,odels 

In the previous section we found that the BEC can have Lagrange invari- 
ants. Let us assume that initially all atoms in the trap are in condensate. 
So, equation (9a) contains the quantity r(r,t) = Up{r,t), which is just the 
rate of pumping, i.e. the rate of increase of the population in ground state 
due to the cooling of atoms from upper-lying mode. Let us find the condi- 
tions for i/p(r, t), under which the pumping of the ground state mode does 
not break the flux lines (i.e. the all local invariants remain to be conserved). 
As it follows from the equation (15), Uc(51 is always a Lagrange invariant 
(see equation (16)) if the gauge condition (14) is satisfied. However, if the 
flux lines are broken so that a vorticity changes, then the condition (14) is 
necessarily broken as well as Uc^l becomes non-invariant. So, we come to 
the simple conclusion: to conserve a given vorticity we must keep the regime 
of pumping to he such one to do the condition (I4) being always valid. Ob- 
viously, the most simple requirement is ^ = a, or as it follows from equation 
(14), ^ 

ix{v,t) = ]^vl{v,t). (17) 

Furthermore, this requirement automatically means that Uc = Vc and the 
invariant v„51 contains the velocity itself. 

For simplicity consider the situation, when the velocity of laser mode 
changes only in given direction. In this case, using equations (9a), (11), and 
(17), we obtain 



^ + i^p{r,t) + V |V2nee„y C/(*™f)(r) +5(ne + 2n) - * = 0, 

_ (18) 

where = Vc/|vc|. Result (18) gives the connection between experimen- 
tally controllable quantities Uc, h, iyp{r,t), and ?7(*™P)(r). Of course, this 
equation should further be solved numerically for particular experimental 
conflgurations. Note that, of course the relation (18) does not solve all 
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problems of atom laser beam stability, but it gives very useful tool for doing 
a choice of the parameters of experimental setup. If the condition (18) is 
satisfied, then at least the problem of angular momentum transfer is solved 
in the frame of the done approximations. 

4-2. Pulsed models 

Here we consider the models of pulsed (not continuously refilled) atom 

laser on the example of rf-transition scheme [15, 16, 40]. In this scheme 
the output coupler includes resonant monochromatic radio frequency field 
transferring atoms in some hyperfine state F from the trapped into un- 
trapped magnetic sublcvcls. In the case of ^^Na atoms F=l, so that s = — 1 
corresponds to the trapped state, s = and s = 1 corresponds to the un- 
trapped and the repelled sublevels, respectively (here F is the total angular 
momentum, s is the magnetic quantum number) . Equation (5) now becomes 
[50] 

= |-|^A + a.a;,/C/i*-^')(r)+5n,(r,i)|l's(r,0 

+nnY.(^s,s'+i + Ss,,'-i)^sir,t), (19) 

s' 

where in rotating wave approximation 

$,(r,t) =e---'-/*(4(r,t)), s,s' e {-1,0, +1}, (20) 

ips{r,t) is the quantum field operator for atoms belong the sublevel with 
given s, cOrf is the frequency of applied resonant rf field, Ug = The 
coupling constant 

1 13 1 

nn = guBohr-^ (21) 

refers to the Rabi frequency due to the rf field. For a small coupling strength 
the process of atoms leaking out of the resonance points is faster than the 
Rabi oscillations. So we further neglect the coupling into state s = +1 since 
it is proportional to O^. Using such approximation and writing the relation 
(7) for each sublevel s we get for the density of atoms in a laser beam no 
the following relation, analogous to the formula (9a): 

^ + V(novo) = 2J^V^sin(^c - ^o), (22) 

where ric = n_i. Equation (22) has a clear physical sense: variation of the 
atom beam density oscillates due to the differences of the condensate and 
beam phases. 
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Equation similar to (9b) looks 

{^Hv„V)}v„.-^. (23, 



where 



Ao = -^^^ + ut'''\r)+gnoir,t) + hn.f^ cosiO, - Oo) (24) 
2m JriQ V n-Q 



is the new chemical potential. From equations (17), (22), (23), and (24) we 
easily obtain the condition similar to (18): 

^ + V I V2noe„ Jt/r'^^)(r) +5no(r,0 - + n^.l^cos{e, - Bo) 

at V 2m A/no V no 



= 2Ov^n^sin(0c - ^o)- (25) 

Note, as it can be seen from the condition (25) the temporal change of the 
beam atoms population depends on the phase difference {9c — Oq), which is 
determined by the frequency LOrf- 

Direct comparison of the formulae (18) and (25) shows that the outcome 
(i.e. laser beam itself) in the continuous case can be stabilized easier than 
in the pulsed regime. 

In conclusion, using the quasi-hydrodynamic approximations we have 
found the local invariants of Bose-Einstein condensate in trapped alkali 
gases. Particularly we obtained the Lagrange invariant, which ensures the 
vorticity to be a freezing-in field. The obtained results can be directly 
applied to the creation of highly coherent atomic beam generators (atom 
lasers) with well controlled angular momentum. Both the pulsed laser and 
laser with continuous refilling are considered. The optimal conditions for 
the pumping modes have been found (again for the both schemes). 
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